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The Evolution of the Visible and Hidden Star Formation 
in the Universe: Implication from the Luminosity 
Functions at FUV and FIR 



T. T. Takeuchi^ V. Buat^, and D. Burgarella^ 

Abstract. Based on GALEX and IRAS/Spitzer datasets, we have found 
that both FUV and FIR luminosity functions (LFs) show a strong evolution 
from z = to 2; = 1, but the FIR LF evolves much stronger than the FUV one. 
Consequently, the FIR/FUV luminosity density ratio increases from 4 (z = 0) 
to 15 {z — 1). It means that more than 80 % of the star-forming activity in the 
Universe is hidden by dust at z = 1. To explore this issue further, we have per- 
formed a combined analysis of the galaxy sample in FUV and FIR. For the Local 
Universe we used GALEX-IRAS sample, whereas at z = 1 we used the Lyman- 
break galaxy sample selected by GALEX bands constructed by Burgarella et 
al. (2005), which is known to be representative of visible (i.e., non-obscured) 
star- forming galaxies at z — 1. From these datasets, we constructed the LFs 
of the FUV-selected galaxies by the survival analysis to, take into account the 
upper-limit data properly. We discovered that the FIR LF of the Lyman-break 
galaxies show a significant evolution comparin g with the loca l FIR LF, but it is 
a factor of 2-3 lower than the global FIR LF (.Le Floc'h et al., 2005 '). This indi- 
cates that the evolution of visible galaxies is not strong enough to explain the 
drastic evolution of the FIR LF. Namely, a FIR-luminous, rapidly diminishing 
population of galaxies is required. 



1. Introduction 

Newly formed massive stars emit strong far-ultraviolet (FUV) radiation. How- 
ever the attenuation of the FUV hght by the interstellar dus t is a major issue 
to derive quantitative SFR from the FUV, even at low-z (e.g.. lBuat et al.ll2005l ) 
On a global point of view, the recent observations conducted by Spitzer and 
GALEX have allowed to build the tot al-IR (TIR) and FUV luminosity func- 
tions and densities from z = to 2: = 1 ("Arno uts et al.ll2005l : ISchiminovich et al 



20051 : iLe Floc'h et al.1 l2005l : iPerez-Gon zalcz eLaLl l2005|) Connecting what is 



seen in FIR and FUV (rest frame) from low- to high- 2: is crucial to understand 
above issues, but a new challenge at the same time. In this work, we try to 
examine which population is responsible for the different evolutions of FUV and 
FIR galaxies through univariate and bivariate luminosity functions (LFs). At 
FUV, we made use of GALEX datasets, and at FIR, we used IRAS data for 
the local sample and Spitzer data for higher-z ones. We use the cosmological 
parameters Hq = 72 km s~^ Mpc~^, Qm = 0.3 and J7a = 0.7. 
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Figure 1. Left: The evolution of the FUV and total-IR (dust) luminosity 
functions obtained from GALEX and IRAS/Spitzer. Right: The evolution 
of the star formation rate density directly obtained from FUV and total IR 
luminosity densities. 



2. Univariate Analysis 

The evolution of the FUV and dust LFs a re shown i n Fig. [H Local LFs at 
FUV and FIR are taken from Wvder et al. I (|2005) and lTakeuchi et al. I ([2003), 



respectively. First, it is worth mentioning the (well-known) difference of the 
local LF shape of FUV and dust: for the dust LF, bright galaxies (L ^ 10^° Lq) 
are much more numerous than those in the FUV. This leads to the difference in 
the main population contributing to the total emitted energy. In the FUV, the 
main contributor is L^, galaxies, with fainter galaxies emitting a non-negligible 
fraction of energy at z > 0.5. In contrast, the effect of the evolution appears 
in the bright end for the dust LF. The contribution from the most luminous 
galaxies increases with redshift. 

By integrating the LFs over luminosities, we obtain luminosity densities at 
FUV and FIR. Both luminosity densities show a significant evolutionary trend, 
but the dust luminosity density evolves much faster than that of the FUV. 
Consequently, the ratio Pdust/PFVY increases toward higher-z, from ~ 4 (local) 
to ~ 15 (z ^ 1), i.e., the dust luminosity dominates the universe at z ~ 1 
(|Takeuchi et al. ll2005al 1. 



We interpret the data in terms of SFR. Ass uming a consta nt SFR over 
10^ yr, and Salpeter initia l mass functioii (IMF ) (jSalpeteil Il955l . mass range: 



0.1-100 Mq), Starburst99 (jLeitherer et al.lll999l ) gives the relation between the 



SFR and L(FUV) = i^L^ at FUV (1530 A), 

log L(FUV) = 9.51 + log SFR . (1) 

For the IR, to transform the dust emission to the SFR, we assume that all the 
stellar light is absorbed by dust. Then, we obtain the following formula under 
the same assumption for both the SFR history and the IMF as those of the 
FUV, 



log L(dust) = 9.75 + log SFR 



(2) 
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However, a significant fraction of the dust emission is due to the heating of grains 
by old stars which is not directly related to the recent SFR. iHira shita et al. 
(|200S ) found that about 30 % of the dust heating in the nearby galaxies comes 
from stars older than 10*^ yr. Adopting this correction, we obtained the evolu- 
tion of the star formation rate densities from FUV and dust (/Osfr(FUV) and 
/OSFR(dust)) which are presented in Fig. [TJ We clearly see that the fraction of 
hidden SFR increases with redshifts, and at z ~ 1 more than 80 % of the cosmic 
SFR is hidden by dust. 



3. Bivariate Analysis: Star-formation LF in the Local Universe 

We could see that most of the star-forming activity in the Universe is hidden by 
dust at z ~ 1. Then, a natural question arises: which population of galaxies is 
due to this evolution? To explore this issue further, we made bivariate LF analy- 
sis at FUV and FIR. First we have carefully constructed FUV and FIR-selected 
samples based on GALEX and IRAS data. Using these datasets, we estimated 
the distribution of the t otal luminosity f rom y oung stars (star-formation lumi- 
nosity). Again following iHirashita et al! ( 20031 ) . we account for the dust heating 



by stars older than 100 Myr. Then the star-formation luminosity is expressed 
as LsF = Lpuv + (1 ~ 'n)^TiR w here r? is the fraction of the TIR emission 
by o ld stars. We adopt rj = 0.3 ( Hirashita et al.l [20031 : llglesias-Paramo et al.l 



20061 ). Therefore LsF 

can be written as Lsf = -^puv + 0.7-Ltir- The star for- 



mation luminosity function from young stars is calculated for each sample using 
the weigh ting method. For details of the upper limit treatment, see 

Buat et al. I (|2006l ). We found that the 60 fim luminosity is a robust tracer of 



the luminosity of young stars, while the FUV flux alone (without any correc- 
tion) misses a large part of the total emission of the FUV selected galaxies. This 
trend is related to the relation found between the lu minosity (or star f ormation 



rate) of the galaxies and their dust attenuation (e.g.. lBuat fc Burgarell a 1998). 
Both luminosity functions are consistent for intermediate luminosities: in the 
nearby universe these galaxies are detected equally well in FIR and in FUV. For 
-^SF > 5 X 10^^ Lq, the star-formation LF issued from the FIR selection is higher 
than that from the FUV one and the discrepancy increases with the luminosity: 
we miss intrinsically bright galaxies which appear much fainter in FUV. Deeper 
wide-area surveys like those of AKARI are waited for further anal ysis. Details 
of the sample construction and analysis are thoroughly explained bv lBuat et al. I 
( 20061 ). Then, how is the situation at higher redshifts? We tried to explore it 



by using a Lyman-break galaxy sample at z ~ 1 in the next section. 



4. Bivariate Analysis: Lyman-break galaxies at z = 1 

Lyman-break galaxies (LBGs) at z ^ 1 are selected by making use of GALEX 
FUV and NUV bands ( Burgarella et al.ll2006l ). It is a suitable method to pick up 
actively star-forming galaxies at these redshifts. To see this, we first estimated 
the LF of LBGs and compared it with the total LF of FUV-selected galaxies at 
z ^ 1 . We have used the L BG sample in the Chandra Deep Field South prepared 
bv [Burgarella et al] ( 20061 ). For the L F estimation, we applied 1/Vmy and C~- 
methods according to the recipes of iTakeuchi. Yoshikawa. Sz Ishiil (j2000l l and 
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Figure 2. The LFs of LBG sample at z = 1. Left: the total FUV LFs of 
FUV-selected galaxies and LBGs. Right: the dust LF of LBGs. Error bars 
represent 2a (95-% CL). 



Takeuchil (|2000l '). The LFs are shown in Fig. El We see that both agree very well 



at the brighter regime, though LBGs are much fewer than purely FUV-selected 
ones at the faint end. 

Then, how much is the contribution of the LBGs to the total dust emission? 
We can address this issue by estimating the dust LF of LBGs. For this, we 
used Spitzer MIPS 24 /im flux to estimate the total dust luminosity. Since at 
z ~ 1 it corresponds to 12 iim, we made use of the conversion formula provided 
by iTakeuchi et al.l ( 2005bl ). The validity of this formu la is guaranteed by a 



comparison with longer wavelength observation at FIR ( Takeuchi et al.ll2006l ) 



In this analysis, we applied the Kaplan-Meier estimator to obtain the dust LF, 
to utilize the information of upper limit data properly. The obtained LF is 
shown in the right panel of Fig. [2j We found that the dust LF of LBGs is a 
factor of 2-3 lower than that of the total dust LF of IR-selected galaxies at 
the same redshifts (jLe Floc'h et al . 2005). This suggests that the evolution of 
visible galaxies is not strong enough to explain the drastic evolution of the FIR 
LF, and a FIR-luminous, rapidly diminishing population of galaxies is required. 
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